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We study inhibitory coherence (i.e., collective coherence by synaptic inhibition) in an ensemble of 
globally-coupled type-I neurons which can fire at arbitrarily low frequencies. No inhibitory coherence 
is observed in a homogeneous ensemble composed of only subthreshold neurons (which cannot fire 
spontaneously without noise). By increasing the fraction of (spontaneously firing) suprathreshold 
neurons P SU pra, heterogeneity- induced inhibitory coherence is investigated in a heterogeneous ensem- 
ble of subthreshold and suprathreshold neurons. As Psupra passes a threshold P 3 * upra , suprathresh- 
old neurons begin to synchronize and play the role of coherent inhibitors for the emergence of 
inhibitory coherence. Thus, regularly-oscillating ensemble-averaged global potential appears for 
Psupra > Psupra- For this coherent case suprathreshold neurons exhibit coherent mixed- mode oscil- 
lations with a fast subthreshold (small-amplitude) hopping frequency and a lower spiking frequency. 
By virtue of their coherent inhibition, sparsely synchronized suprathreshold neurons suppress noisy 
activities of subthreshold neurons. Thus, only coherent subthreshold hoppings appear in the indi- 
vidual potentials of subthreshold neurons. We also characterize the inhibitory coherence in terms 
of the "statistical-mechanical" spike-based and correlation-based measures and find that the de- 

> Psupra- Finally, effect 



gree of inhibitory coherence increases with increasing Psupra for Psupra 
of sparse randomness of synaptic connectivity on the inhibitory coherence and universality of the 
heterogeneity-induced inhibitory coherence are briefly discussed. 

PACS numbers: 87.19.1m, 87.19.1c 



I. INTRODUCTION 



Recently, much attention has been paid to rhythms of 
the brain Coherence of neural oscillations may be 
used for efficient sensory and cognitive processing (e.g., 
feature integration, selective attention, working memory, 
and decision making) 0, H| . This kind of neural coher- 
ence is also correlated with pathological rhythms asso- 
ciated with neural diseases (e.g., epileptic seizures and 
tremors in the Parkinson's disease) Here, we are in- 
terested in these coherent brain rhythms. A brain circuit 
is composed of a few types of excitatory principal cells 
and diverse types of inhibitory interneurons. Interneu- 
ron diversity increases the computational power of prin- 
cipal cells [lj. Effect of chemical synapses on coherent 
brain rhythms has been much investigated in neural sys- 
tems composed of excitatory and/or inhibitory neurons 
0, @- Historically, recurrent excitation between prin- 
cipal cells is the conventional coherence mechanism [6|. 
However, when the decay time of the synaptic interaction 
is enough long, mutual inhibition between interneurons 
(rather than excitation) may synchronize individual neu- 
ral firings [7, 8]. By providing a coherent oscillatory out- 
put to the principal cells, interneuronal networks play the 
role of the backbones (i.e., pacemakers) of many brain 
rhythms such as the thalamocortical spindle rhythms 
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[j| [l(| and the fast gamma rhythms in the hippocampus 
and the neocortex |ll| - |l4| . When the feedback between 
the excitatory and the inhibitory populations is strong, 
neural coherence occurs via the "cross-talk" between the 
two populations [l3rll6| . In these computational stud- 
ies of neural coherence, different types of network ar- 
chitectures have been considered Q; all-to-all networks 
where every neuron is coupled to every other neuron, 
sparse random networks where synaptic connections are 
sparse, and complex networks such as small- world net- 
works (with predominantly local connections and rare 
long-distance connections) [ijj and scale-free networks 
(with a few percent of hub neurons with an exceptionally 
large number of connections) [TH . 

Neurons in the neural system exhibit a variety of mor- 
phological and physiological properties. However, close 
to threshold, this remarkable richness may be grouped 
broadly into two basic types of excitability, often referred 
to as type I and type II [19j . When the strength of a con- 
stant input current passes a threshold, type-I neurons can 
fire at arbitrarily low frequencies and they can smoothly 
encode the strength of the input into the output firing 
frequency. In contrast, type-II neurons have a non-zero 
minimum frequency of firing and they fire in a narrow 
frequency band which is relatively insensitive to changes 
in the strength of the applied current. Different types of 
excitability occur because neurons have different bifur- 
cations of resting and spiking states 20]. For the type 
I neurons, oscillations emerge via a saddle-node bifurca- 
tion on an invariant circle. As the bifurcation parameter 
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(i.e., strength of the injected current) passes a threshold, 
the stable and the unstable fixed points coalesce and then 
disappear, leaving a large-amplitude stable periodic or- 
bit. This is a global bifurcation and the frequency of the 
global loop can be arbitrarily small. On the other hand, 
for type-II neurons a transition from a resting state to a 
periodically spiking state occurs through a Hopf bifurca- 
tions with a finite non-zero firing frequency. According 
to their bifurcations, neurons may also be classified into 
integrators and resonators [2l|. Type-I neurons act as 
integrators without subthreshold oscillations, and they 
prefer high-frequency input: the higher the frequency of 
the input, the sooner they fire. In contrast, type-II neu- 
rons exhibit damped subthreshold oscillations and act as 
resonators: they prefer oscillatory input with the same 
frequency as that of damped oscillations. According to 
their excitability type, neurons make distinctly different 
responses to stimuli which have important implications 
for their distinct roles in generating population rhythms 
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In this paper, we study inhibitory coherence (i.e., col- 
lective coherence by synaptic inhibition) in an ensemble 
of globally-coupled type-I neurons. Neural models ex- 
hibiting the type-1 excitability include the Connor model 
for the crab leg axons 12711 . the Wang-Buzsaki model for 
inhibitory interneurons the Hindmarsh-Rose model 
[28| . and the Morris-Lecar (ML) model [29| under some 
circumstances. In Section UH we describe the biologi- 
cal conductance-based ML neuron model with voltage- 
gated ion channels. The ML neurons (used in our study) 
exhibit the type-I excitability, and they interact via in- 
hibitory GABAergic synapses whose activity increases 
fast and decays slowly. Inhibitory coherence (which is our 
main concern) is particularly important because it plays 
a significant role in integration of sensory and cognitive 
information; for example, impaired inhibitory coherence 
is believed to be associated with schizophrenia and at- 
tention deficit disorder [3(3, Hl|. Hence, it is important to 
understand mechanisms for the emergence of inhibitory 
coherence. Many works exploring mechanisms of neu- 
ral coherence were done in neural systems composed of 
spontaneously firing (i.e., self-oscillating) suprathreshold 
neurons (above the threshold) 0, Q • For this case, neu- 
ral coherence occurs via cooperation of regular firings of 
suprathreshold neurons. Unlike the suprathreshold case, 
subthreshold neurons (below the threshold) cannot fire 
spontaneously without noise; they can fire only with the 
help of noise. Stochastic excitatory coherence (i.e., collec- 
tive coherence between noise-induced spikings by synap- 
tic excitation) was observed in a population of excita- 
tory subthreshold neurons [H, [33| . Due to the stochas- 
tic excitatory coherence, synaptic current, injected into 
each individual neuron, becomes temporally coherent. 
Hence, temporal coherence resonance of an individual 
subthreshold neuron in the network may be enhanced. 
Furthermore, stochastic inhibitory coherence (i.e., collec- 
tive coherence between noise-induced spikings by synap- 
tic inhibition) was also investigated in a population of 



inhibitory subthreshold ML neurons exhibiting the type- 
II excitability [34j . Weak stochastic inhibitory coherence 
was thus found to appear via cooperation of individ- 
ual irregular oscillations (i.e., a regular small-amplitude 
ensemble-averaged oscillation emerges from sparsely syn- 
chronized neurons discharging irregularly at lower rates 
than the network oscillation). These sparsely synchro- 
nized neural oscillations have been intensively investi- 
gated in other types of neural networks 35] and they 



are believed to be associated with cortical rhythms in 
cognition [e.g. ultrafast rhythm (100-200 Hz), gamma 
rhythm (30-100 Hz) and beta rhythm (15-30 Hz)] with 
irregular and sparse neural discharges 0, |35[ . 

In contrast to the case of subthreshold type-II ML 
neurons, no stochastic inhibitory coherence is observed 
in a homogeneous population of subthreshold type-I 
ML neurons. Hence, subthreshold type-I integrator 
neurons (without subthreshold oscillations) seem to be 
much more difficult to synchronize by inhibition than 
subthreshold type-II resonator neurons (exhibiting sub- 
threshold oscillations). To take into consideration the 
effect of (spontaneously firing) suprathreshold neurons 
on the inhibitory coherence, we consider a heterogeneous 
inhibitory ensemble of subthreshold and suprathreshold 
type-I ML neurons. Heterogeneity (or diversity) has been 
found to make constructive effects on collective coher- 
ence in various physical, biological, neural, and social 
systems (36|. In Section Hill we investigate heterogeneity- 
induced inhibitory coherence by increasing the fraction of 
suprathreshold neurons P SU pra in the whole population. 
As Psupra passes a threshold value Ps* upra suprathresh- 
old neurons begin to synchronize and they play the role 
of coherent inhibitors for the emergence of inhibitory co- 
herence in the whole heterogeneous population. Thus, 
for Psupra > Psupra the ensemble-averaged global po- 
tential Vq exhibits a regular small-amplitude oscillation. 
For this coherent case, individual suprathreshold neurons 
exhibit intermittent spikings phase-locked to Vg at ran- 
dom multiples of the period of Vq- Due to the stochas- 
tic spike skipping of suprathreshold neurons, the inter- 
spike interval (ISI) histogram has multiple peaks and 
partial occupation occurs in the raster plot of neural 
spikes. In addition to the coherent intermittent spik- 
ings, coherent subthreshold (small-amplitude) hopping 
oscillations also appear in the individual potentials of 
suprathreshold neurons. Thus, sparsely synchronized 
suprathreshold neurons exhibit coherent mixed-mode os- 
cillations with two well-separated frequencies, a fast sub- 
threshold hopping frequency imposed by the collective 
network oscillation and a lower firing frequency of indi- 
vidual suprathreshold neurons. By virtue of their co- 
herent inhibition sparsely synchronized suprathreshold 
neurons suppress noisy activities of subthreshold neu- 
rons. Thus, only coherent fast subthreshold hopping 
oscillations (without spikings) appear in the individual 
potentials of subthreshold neurons. We also characterize 
this heterogeneity-induced inhibitory coherence in terms 
of "statistical-mechanical" spike-based and correlation- 
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based measures, and find that the degree of inhibitory 
coherence increases as P BU pra is increased from Ps Upra - 
fn a real brain, each neuron is coupled to only a certain 
number of neurons which is much smaller than the total 
number of neurons. The effect of sparseness of synaptic 
connectivity on the inhibitory coherence is briefly dis- 
cussed by varying the average number of synaptic in- 
puts per neuron M syn in a random network. Emer- 
gence of inhibitory coherence is thus found to persist un- 
til M syn is larger than a threshold value M* yn . We also 
confirm the universality of the heterogeneity-induced in- 
hibitory coherence in a population of canonical type-I 
quadratic integrate-and-fire neurons 37|. This kind of 
heterogeneity-induced weak inhibitory coherence might 
be associated with cortical rhythms with stochastic and 
sparse neural discharges which contribute to cognitive 
functions in the cerebral cortex (e.g., information inte- 
gration, working memory, and selective attention) [3,l35j. 
Finally, a summary is given in Section [IVI 



II. HETEROGENEOUS ENSEMBLE OF 
INHIBITORY TYPE-I ML NEURONS 

In this section we describe the biological neuron model 
used in our computational study. We consider a hetero- 
geneous inhibitory ensemble of N globally-coupled sub- 
threshold and suprathreshold type-I neurons. As an ele- 
ment in our neural system, we choose the conductance- 
based ML neuron model, originally proposed to describe 
the time-evolution pattern of the membrane potential for 
the giant muscle fibers of barnacles [2^] . The population 
dynamics in this neural network is governed by the fol- 
lowing set of differential equations: 
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Iion,i — Ica.i + Ik,% + Ih.i (2a) 
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mao (v) = 0.5[l + tanh{(v-Vi)/y 2 }], (2d) 

w^v) = 0.5[l+tanh{(w- V 3 )/V 4 }}, (2e) 

r R (v) = l/cosh{>- V 3 )/(2Vi)}, (2f) 

sM = l/[l + e-^- v ^ 5 }. (2g) 

Here, the state of the ith neuron at a time t (measured in 
units of ms) is characterized by three state variables: the 



membrane potential V{ (measured in units of mV), the 
slow recovery variable Wi representing the activation of 
the K + current (i.e., the fraction of open K + channels), 
and the synaptic gate variable Si denoting the fraction 
of open synaptic ion channels. In Eq. (jTJa,) , C represents 
the capacitance of the membrane of each neuron, and the 
time evolution of u, is governed by four kinds of source 
currents. 

The total ionic current Ii on ,% of the ith neuron con- 
sists of the calcium current Ica,ii the potassium current 
Ik,Ii and the leakage current Il,%- Each ionic current 
obeys Ohm's law. The constants gc a , 9k, and are 
the maximum conductances for the ion and the leakage 
channels, and the constants Vc a , Vk, and Vl are the 
reversal potentials at which each current is balanced by 
the ionic concentration difference across the membrane. 
Since the calcium current Ica.i changes much faster than 
the potassium current Ik,u the gate variable rrn for the 
Ca 2+ channel is assumed to always take its saturation 
value rriooivi). On the other hand, the activation variable 
Wi for the K + channel approaches its saturation value 
Wca{vi) with a relaxation time TR{vi)/4>, where tr has a 
dimension of ms and <p is a (dimensionless) temperature- 
like time scale factor. 

Each ML neuron is also stimulated by a DC cur- 
rent Idc,i and a Gaussian white noise £j [see the 2nd 
and 3rd terms in Eq. (JTJi)] satisfying = and 

— $ij <K* — *')> wriere (• • • ) denotes the en- 
semble average. The noise £j is a parametric one which 
randomly perturbs the strength of the applied current 
Idc,ii and its intensity is controlled by the parame- 
ter D. Depending on the system parameters, the ML 
neuron may exhibit either type-I or type-II excitabil- 
ity [2^| . Throughout this paper, we consider the case 
of type-I excitability where gc a = 4 mS/cm 2 , gx — 
8 mS/cm 2 , g L = 2 mS/cm 2 , V Ca = 120 mV, V K = 
-84 mV, V L = -60 mV, C = 20 /iF/cm 2 , cj> = 1/15, 
Vi = -1.2 mV, V-2 = 18 mV, V 3 = 12 mV, and V 4 = 
17.4 mV. (For comparison, a result on the order pa- 
rameter is given for the type-II case where the values 
of the above parameters are the same as those in the 
type-I case except that gc a — 4.4 mS/cm 2 , cfr — 0.04, 
V3 = 2 mV, and V4 = 30 mV.) For the type-I case, a 
transition from a resting state to a spiking state occurs 
for l* DC = 40 [lA/cm 2 via a saddle- node bifurcation on 
an invariant circle [2fJ, and a firing begins at arbitrar- 
ily low frequency. On the other hand, a type-II neuron 
exhibits a jump from a resting state to a spiking state 
through a subcritical Hopf bifurcation for 1* DC h = 93.9 
/iA/cm 2 by absorbing an unstable limit cycle born via 
fold limit cycle bifurcation for I* DC ; = 88.3 /jA/cm 2 [20j | . 
and hence the firing frequency begins from a non-zero 
value. Here a spread in the value of the DC input cur- 
rent I DC is taken into consideration, and thus for each 
subthreshold (suprathreshold) neuron a value of Idca is 
randomly chosen with a uniform probability in the range 

of ( r DC- A > r Dc) i( I DC> I hc + A )} for the type" 1 case and 
in the range of (/* c , - A, J* c ,) \(I* DC h , I* DC ,h + A )] for 
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the type-II case, where the value of the spread parameter 
A is set as A = 10 ^A/cm 2 . 

We consider a heterogeneous inhibitory ensemble of 
N globally-coupled subthreshold and suprathreshold ML 
neurons where the fraction of suprathreshold neurons 
is given by P supra = Ns ^ ra (N supra : number of 
suprathreshold neurons) . The last term in Eq. (JTJl) repre- 
sents the synaptic coupling between neurons in the net- 
work. Each neuron is connected to all the other ones 
through global synaptic couplings. I syn .i of Eq. ([5b) rep- 
resents such synaptic current injected into the ith neuron. 
Here the coupling strength is controlled by the parame- 
ter J and V syn is the synaptic reversal potential. We use 
V syn = —80 mV for the inhibitory synapse. The synaptic 
gat e variable s obeys the 1st order kinetics of Eq. ((Tb) 
fiol [ll| . Here the normalized concentration of synaptic 
transmitters Soc(v), activating the synapse, is assumed 
to be an instantaneous sigmoidal function of the mem- 
brane potential with a threshold t>t in Eq. ((2g) , where we 
set vt = mV and 6 = 2 mV. The transmitter release 
occurs when the neuron emits a spike (i.e., its potential v 
is larger than v t ). For the inhibitory GABAergic synapse 
(involving the GABAa receptors), the synaptic channel 
opening rate, corresponding to the inverse of the synap- 
tic rise time r r , is a = 10 ms -1 , and the synaptic closing 
rate ft, which is the inverse of the synaptic decay time Td, 
is (3 = 0.1 ms -1 [l6j]. Hence, I syn rises fast and decays 
slowly. 

Numerical integration of Eq. (fTJ) is done using the Hcun 
method [38| (with the time step At = 0.01 ms) similar 
to the second-order Runge-Kutta method, and data for 
(Vi, Wi, Si) (i = 1, . . . , N) are obtained with the sampling 
time interval At = 1 ms. For each realization of the 
stochastic process in Eq. (TTJ), we choose a random ini- 
tial point [vi(0),Wi(0),Si(0)] for the ith (i = 1,...,N) 
neuron with uniform probability in the range of w, (0) 6 
(-70,50), u)i(0) G (0.0,0.6), and s 4 (0) e (0.0,1.0). 



III. HETEROGENEITY-INDUCED 
INHIBITORY COHERENCE 

In this section, we are concerned about inhibitory 
coherence in a heterogeneous ensemble of N globally- 
coupled subthreshold and suprathreshold type-I ML neu- 
rons. By increasing the fraction of suprathreshold neu- 
rons P SU pra , we investigate the heterogeneity-induced in- 
hibitory coherence. (Hereafter, for convenience we omit 
the dimensions of Idc, D, and J.) 

We first consider a homogeneous population (corre- 
sponding to the case of P SU pra = 0) composed of only sub- 
threshold type-I ML neurons, and study the inhibitory 
coherence by varying both the coupling strength J and 
the noise intensity D. Emergence of inhibitory coherence 
may be well described by the (ensemble-averaged) global 




FIG. 1: Plots of the order parameter O versus (a) both the 
coupling strength J and the noise intensity D and versus (b) 
D for J = 20 mS/cm in N globally-coupled inhibitory sub- 
threshold type-I ML neurons. The value of Idc,% for each 
subthreshold type-I neuron is randomly chosen with a uniform 
probability in the range of (loc ~ A, Ihc) where Ip C = 40 
/iA/cm 2 and A = 10 nA/cm 2 . (c) Plots of O versus D for 
J = 3 mS/cm 2 in iV globally coupled inhibitory subthreshold 
type-II ML neurons. The value of Ioc,i for each subthreshold 
type-II neuron is randomly chosen with a uniform probability 
in the range of (Inai ~ A, Idc.i) where Ino,i = 88-3 /iA/cm 2 
and A = 10 /iA/cm 2 . 

potential, 

1 N 

In the thermodynamic limit (N — > oo), a collective state 
becomes coherent if AVc(t) (= Vc{t) — Vc(t)) is non- 
stationary (i.e., an oscillating global potential Vg appears 
for a coherent case), where the overbar represents the 
time average. Otherwise (i.e., when AVq is stationary), 
it becomes incoherent. Thus, the mean square deviation 
of the global potential Vq (i.e., time-averaged fluctua- 
tions of Vg), 



O = (V G (t) - V G (t))\ (4) 

plays the role of an order parameter used for describ- 
ing the coherence-incoherence transition [39]. For the 
coherent (incoherent) state, the order parameter O ap- 
proaches a non-zero (zero) limit value as N goes to the 
infinity. Figure Ufa) shows plots of the order parameter 
versus both the coupling strength J and the noise inten- 
sity D. As N is increased, the order parameter tends to 
decrease, independently of J and D. An example of the 
order parameter is shown in Fig.QJb) for J = 20. For any 
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FIG. 2: Order parameters, raster plots of neural spikes, and 
time series of global potentials in the heterogeneous ensem- 
ble of N globally-coupled inhibitory type-I ML neurons for 
J = 20 mS/cm 2 and D = S fiA ■ ms 1/2 /cm 2 ; N = 10 3 in 
(bl)-(b5) and (cl)-(c3). The value of Inc,i for each subthresh- 
old (suprathreshold) type-I ML neuron is randomly chosen 
with a uniform probability in the range of (Ioc ~ A, Ihc) 
[(Ihc^hc + A)] where I* DC = 40 fiA/cm 2 and A = 10 
/iA/cm 2 . Plots of the order parameter O versus the fraction 
of suprathreshold neurons P S u P ra in (al) the whole population 
and in the two subpopulations of (a2) the suprathreshold and 
(a3) the subthreshold neurons. Raster plots and time series of 
the global potential Va in the whole population for Psupra = 
(bl) 0, (b2) 0.2, (b3) 0.4, (b4) 0.6, and (b5) 1.0. Time series 
of the subensemble-averaged potentials V SU pra and V su b in the 
two subpopulations of the suprathreshold and the subthresh- 
old neurons for P 3up ra = (cl) 0.2, (c2) 0.4, and (x c3) 0.6. 
Vertical dashed lines in (cl)-(c3) represent the times at which 
local minima of Vo appear. 



given D, O is found to decrease as N is increased. Hence, 
only incoherent states exist, irrespectively of D. This is 
in contrast to the case of subthreshold type-II neurons 
exhibiting inhibitory coherence. Figure [He) shows plots 
of the order parameter versus the noise intensity for the 
type-II case of J = 3. Unlike the type-I case coherent 
states exist in an intermediate range of noise intensity 
[Df(~ 10.3) < D < D* h {~ 27.9)] where the order param- 
eter approaches a non-zero limit value as N increases. 

To take into consideration the effect of (spontaneously 
firing) suprathreshold neurons on the inhibitory coher- 
ence, we consider a heterogeneous population consist- 
ing of subthreshold and suprathreshold type-I ML neu- 
rons for J = 20. For convenience, we set the value of 
noise intensity as D = 8 and investigate the heterogene- 
ity effect on the inhibitory coherence by increasing the 
fraction of suprathreshold neurons P SU pra- Figure [Hal) 



shows plots of the order parameter O versus P sup ra m 
the whole population. As P sup ra passes a threshold value 
Psuprai— 0-16), a transition from an incoherent to a co- 
herent state occurs. As shown in Fig. [2Jal), it is enough 
to consider only the case of N = 10 3 for the study of 
inhibitory coherence because the order parameter O be- 
comes saturated for N — 10 3 . Hereafter we set the num- 
ber of neurons as N = 10 3 in all cases except the cal- 
culation of the order parameter. For an incoherent case 
of Psupra — 0, the raster plot consists of randomly scat- 
tered sparse spikes and the global potential Vq exhibits 
nearly a stationary irregular oscillation [see Fig. [2fbl)]; 
the amplitude of Vq decreases with further increase in N. 
However, when passing the threshold P t 



supra 



partially- 
occupied "stripes" (composed of spikes and indicating 
collective coherence) appear in the raster plot together 
with regularly-oscillating small-amplitude Vq with fre- 
quency fc (= 13.8 Hz) [see Fig. [2^b2)]. As P sup ra is 
further increased, both the pacing degree of spikes and 
the amplitude of Vq (representing the degree of collective 
coherence) increase, as shown in Figs. [SJb3)-[2Ib5) where 
f G = (b3) 14.3 Hz, (b4) 13.6 Hz, and (b5) 14.2 Hz. This 
kind of weak inhibitory coherence also occurs in each sub- 
population of the subthreshold and the suprathreshold 
neurons. As in the case of the whole population, emer- 
gence of inhibitory coherence in the subpopulations may 
be well described by the subensemble-averaged potentials 

Vsupra and Vsubi 



v.. 



■pro (0 — 



supra i—\ 



V sub {t) 



N. 



sub 



J2 



(5a) 



(5b) 



where N supra (N suo ) is the number of suprathreshold 
(subthreshold) neurons. Then the order parameters 
O supra and Osub, defined by the mean square deviation 

of Vsupra 

and Vsub, 



O. 



supra 



(Vsupra Vsupra 



{Vsub{t)-Vsub{t)f 



(6a) 
(6b) 



may be used for describing the coherence-incoherence 
transitions in the subpopulations of suprathreshold and 
subthreshold neurons, respectively. Plots of O SU pra and 
Osub versus P SU pra are shown in Figs. [2ja2) and[2ja3), 
respectively. Coherent transition in each subpopulation 
occurs at the same threshold value Pguprai— 0.16). For 
the case of coherent states, not only V SU pra but also V su b 
exhibits regular oscillations whose amplitudes increase 
as Psupra is increased [see Figs. EJcl)-[2]Jc3) where verti- 
cal dashed lines in V SU pra and V su b denote the times at 
which local minima of Vq appear] . Both V SU pra and V su b 
are phase-locked to Vq. 

To further understand the emergence of inhibitory co- 
herence, we examine the individual and the global out- 
put signals in the subpopulations of subthreshold and 
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FIG. 3: Time series of the individual and the global poten- 
tials, the average firing probability, and the interspike interval 
(ISI) histogram in the heterogeneous ensemble of N(= 10 3 ) 
globally-coupled inhibitory type-I ML neurons for J = 20 
mS/cm 2 and D = 8 /j,A ■ ms 1//2 /cm 2 . The value of Iuc,% 
for each subthreshold (suprathreshold) type-I ML neuron is 
randomly chosen with a uniform probability in the range of 
(P DC - A, r DC ) [{IocJ'dc + A)] where P DC = 40 M/cm 2 
and A = 10 /xA/cm 2 . The individual potential v\ of the first 
neuron and the global potential Vsupra in the subpopulation 
of suprathreshold neurons for P SU pra = (al) 0.1, (a2) 0.2, (a3) 
0.4, (a4) 0.6, and (a5) 1.0. Vertical dashed lines in (a2)-(a5) 
represent the times at which local minima of Vsupra appear, 
(b) Plot of the average firing probability Pf tSU b versus the frac- 
tion of suprathreshold neurons P supra . The individual poten- 
tial vi of the first neuron and the global potential V su b in the 
subpopulation of subthreshold neurons for Psupra = (cl) 0, 
(c2) 0.1, (c3) 0.2, (c4) 0.4, and (c5) 0.6. Vertical dashed lines 
in (c3)-(c5) denote the times at which local minima of V su b 
appear. ISI histograms in the whole population for P 3U pra 
(dl) 0, (d2) 0.2, (d3) 0.4, (d4) 0.6, and (d5) 1.0; each ISI 
histogram is composed of 5 x 10 4 ISIs and the bin size for the 
histogram is 5 ms. Vertical dotted lines in (d2)-(d5) denote 
integer multiples of Tg (period of Vg)- 



suprathreshold neurons. Figures EJal )-[3]^a5) show the 
time series of the individual potential v\ of the first neu- 



ron and the time series of the global potential V s 



supra 



the subpopulation of the suprathreshold neurons. V s 



supra 



exhibits a regular small-amplitude oscillation for a co- 
herent case, while it shows a nearly stationary irregu- 
lar oscillation for an incoherent case. For the case of 
coherent states, individual suprathreshold neurons ex- 



hibit intermittent spikings phase-locked to V supra at ran- 
dom multiples of the period of V supra [see Figs. [3{a2)- 
[3ja5) where dashed lines denote the times at which local 
minima of V supr a appear]. This "stochastic phase lock- 
ing" leading to stochastic spike skipping is well shown in 
the ISI histogram with multiple peaks [see Figs. [3jcl2)- 
[3£d5)]. In addition to these coherent intermittent spiking 
phases, coherent subthreshold (small-amplitude) hopping 
oscillations also appear in the individual potentials of 
suprathreshold neurons. Thus, suprathreshold neurons 
exhibit coherent mixed-mode oscillations with two well- 
separated frequency scales, a fast subthreshold hopping 
frequency f n imposed by the collective network oscilla- 
tion with frequency /g (— 14 Hz) and a lower spiking 
frequency f s of individual suprathreshold neurons; f s = 
(a2) 3.6 Hz, (a3) 2.8 Hz, (a4) 2.3 Hz, and (a5) 1.8 Hz. 
These sparsely synchronized suprathreshold neurons play 
the role of coherent inhibitors for the emergence of in- 
hibitory coherence in the whole heterogeneous popula- 
tion, as shown below. On the other hand, for the in- 
coherent case only stochastic intermittent spikings occur 
without any coherent hoppings, as shown in Fig. [3jal). 
Figure [2Jb) shows the plot of the average firing proba- 
bility Pf,sub versus P supra in the subpopulation of sub- 
threshold neurons (i.e., time-averaged fraction of firing 
subthreshold neurons in the subpopulation of subthresh- 
old neurons). Due to inhibition Pf jSU b decreases dramat- 
ically with respect to P supra - For P SU pra > 0.02, one 
can disregard spikings of subthreshold neurons because 
Pf,sub{~ 10 -8 ) becomes very small. The time series of 
the individual potential v-y of the first neuron and the 
time series of the global potential V su b m the subpopula- 
tion of the subthreshold neurons are shown in Fig. [Ucl)- 
E£c5). A regular oscillation with small amplitude oc- 
curs in V su b for a coherent case, while V su b exhibits a 
nearly stationary irregular oscillation for an incoherent 
case. For the case of coherent states, sparsely synchro- 
nized suprathreshold neurons suppress noisy activities of 
subthreshold neurons by virtue of their coherent inhibi- 
tion, and then individual subthreshold neurons exhibit 
only coherent subthreshold hoppings (without spikings), 
in contrast to the suprathreshold case. Figures Eljdl)- 
[3fd5) show the ISI histograms in the whole population. 
(As shown above, spiking neurons in the whole popu- 
lation are just suprathreshold ones for P SU pra > 0.02.) 
The ISI histogram for P swpra = has a very long tail, 
and hence the average value (ISI)(~ 4926 ms) of ISIs 
is very large. As Psupra passes the threshold f s * upr . a , 
multiple peaks tend to appear at integer multiples of 
T G (period of Vg) [i.e., nT G (n = 1,2,3,...)] [e.g., see 
the ISI histogram for P supr a = 0.2; vertical dotted lines 
in the histogram denote integer multiples of Tq{= 72 .4 
ms)]. As Psupra is further increased, ISI histograms with 
more distinct multiple peaks appear due to the stochastic 
spike skipping of the suprathreshold neurons, as shown 
in Figs.[3];d3)-[3i;d5) where T G = (d3) 69.9 ms, (d4) 73.3 
ms, and (d5) 70.3 ms. The most probable peak appears 
at 2 Tq, and hence suprathreshold neurons fire mostly in 
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alternate global cycles. 

We characterize the heterogeneity-induced inhibitory 
coherence in terms of two kinds of "statistical- 
mechanical" spike-based and correlation-based measures. 
As shown in Figs. [2^bl)-[2][b5), inhibitory coherence may 
be well visualized in the raster plot of spikes. For a coher- 
ent case, the raster plot is composed of partially-occupied 
stripes (indicating collective coherence). To measure the 
degree of the collective coherence seen in the raster plot, 
a new spike-based measure M s was introduced by consid- 
ering the occupation pattern and the pacing pattern of 
neural spikes in the "stripes" [34| . Particularly, the pac- 
ing degree between spikes is determined in a statistical- 
mechanical way by quantifying the average contribution 
of microscopic individual spikes to the global potential 
Vg- The spiking coherence measure M, of the ith stripe 
is defined by the product of the occupation degree Oi of 
spikes (representing the density of the ith stripe) and the 
pacing degree Pi of spikes (denoting the smearing of the 
ith stripe): 



Mi=Oi-Pi. 



(7) 



The occupation degree Oi in the ith stripe is given by 
the fraction of spiking neurons: 



Oi = 



N 



0) 



N 



(8) 
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FIG. 4: "Statistical-mechanical" coherence measures in the 
heterogeneous ensemble of N(= 10 3 ) globally-coupled in- 
hibitory type-I ML neurons for J = 20 mS/cm 2 and D = 8 
fiA ■ ms ly/2 /cm 2 . (a) Spike-based coherence measure M a : (al) 
plot of the average occupation degree (Oi) versus the fraction 
of suprathreshold neurons P 3upra , (a2) plot of the average 
pacing degree (Pi) versus P S u P ra, and (a3) plot of the spiking 
coherence measure M s versus P aup ra- To obtain (Oi), (Pi), 
and M s , we follow the 3 x 10 stripes for each P S u pra - (b) 
Correlation-based coherence measure M c : (bl) plot of M c 
versus P S u pra in the whole population and plots of Mc Supra ^ 
and Mc aub ^ versus P S u P ra in the two subpopulations of (b2) 
suprathreshold and (b3) subthreshold neurons. The number 
of data used for the calculation of a cross-correlation function 
for each P S u P ra is 2 12 . 



(s) 

where N£ is the number of spiking neurons in the ith 
stripe. For the full occupation, Oi = 1, while for the 
partial occupation Oi < 1. The pacing degree Pi of each 
microscopic spike in the ith stripe can be determined in 
a statistical-mechanical way by taking into consideration 
its contribution to the macroscopic global potential Vg- 
Each global cycle of Vg begins from its left minimum, 
passes the central maximum, and ends at the right mini- 
mum; the central maxima coincide with centers of stripes 
in the raster plot [see Figs. HJb2)-[2Ib5)]. An instanta- 
neous global phase <J>(i) of Vg is introduced via linear 
interpolation in the two successive subregions forming 
a global cycle. The global phase $(i) between the left 
minimum (corresponding to the beginning point of the 
ith global cycle) and the central maximum is given by: 



= 2n{i - 3/2) + 7T 



t - t 



(mm) 



.(max) Amin) J 
si ~ l i J 

for 4 mm) < t < t[ max) (i = 1, 2, 3, . . . ), (9) 

and $(t) between the central maximum and the right 
minimum (corresponding to the beginning point of the 
(i + l)th cycle) is given by: 



$(t) = 2tt(z - 1)+tt 



t-t 



(max) 



Amin) .(max) 



where ^ mm ) j s the beginning time of the ith global cycle 
(i.e., the time at which the left minimum of Vg appears in 
the ith global cycle) and ^■■ max ^ is the time at which the 
maximum of Vg appears in the ith global cycle. Then, 
the contribution of the kth microscopic spike in the ith 

( s) 

stripe occurring at the time V k ' to Vg is given by cos , 
where <£>/,. is the global phase at the fcth spiking time 
[i.e., $fc = 3>(tj?')]. A microscopic spike makes the most 
constructive (in-phase) contribution to Vg when the cor- 
responding global phase is 27m (n = 0, 1, 2, . . . ), while 
it makes the most destructive (anti-phase) contribution 
to Vg when $i is 2n(n — 1/2). By averaging the contri- 
butions of all microscopic spikes in the ith stripe to Vg, 
we obtain the pacing degree of spikes in the ith stripe, 



P 



1 Sl 

-^COS$fe, 

fe=l 



S, 



(11) 



where Si is the total number of microscopic spikes in the 
ith stripe. By averaging M, of Eq. (J7J) over a sufficiently 
large number N s of stripes, we obtain the spike-based 
coherence measure AL: 



1 Ns 
M s = — > A'U 



(12) 



for t'f nax ' > <t< t^™" 1 ^ (i = 1, 2, 3, . . . ), (10) By varying P sup ra, we follow 3 x 10 3 stripes and mea- 
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sure the degree of collective spiking coherence in terms 
of (Oi) (average occupation degree), (P,-) (average pac- 
ing degree), and M s for 13 values of P sup ra in the co- 
herent regime, and the results are shown in Figs. EJal)- 
[4](a3). As P sup ra is increased, the average occupation de- 
gree {Oi) (denoting the average density of stripes in the 
raster plot) increases slowly, but its values are very small 
((Oi) < 0.05); only a fraction (less than 1/20) of the total 
neurons fire in each stripe [see Figs. [2jb2)-[2][b5)]. This 
partial occupation results from stochastic spike skipping 
of individual neurons seen well in the multi-peaked ISI 
histograms [see Figs. GUd2)-[3Jd5)]. On the other hand, 
the average pacing degree (Pi) increases rapidly near the 
threshold P£ upra , and then it grows slowly. This tendency 
may be understood from the change in the structure of where 
the ISI histograms. As P SU pra is increased, clear well- 
separated multiple peaks appear, and hence the average 
pacing degree of the stripes becomes better with increas- 
ing Psupra- In most region of the coherent region, the 
values of (Pi) are large in contrast to (Oi). However, the 
spiking measure M s of Eq. (|12p (representing the col- 
lective spiking coherence) is very low due to the partial 
occupation in the raster plot. 

For the coherent case, subthreshold neurons exhibit 
only coherent subthreshold hoppings without spikings. 
Hence, only suprathreshold neurons (exhibiting coherent 
intermittent spikings) make contribution to the spike- 
based measure M s . From now on, we use another 
statistical-mechanical measure based on the ensemble av- 
erage of cross-correlations between the global potential 
and the individual potentials [4(j> and measure the de- 
gree of the inhibitory coherence in the subpopulations 
of the subthreshold and the suprathreshold neurons as 
well as in the whole population. The inhibitory coher- 
ence in the whole population is quantified in terms of the 
coherence measure M c given by the ensemble average of 
the global- individual cross-correlations Ci (0) between V G 
and Vi at the zero-time lag: 



M r 



1 N 



(13) 



i=l 



Here the normalized cross-correlation function Ci (r) be- 
tween V G and Vi is given by 



Ci(r) 



AV G (t + r) Avi(t) 



(14) 



where r is the time lag, AV G (t) = V G (t) - V G (t), 
Avi(t) — Vi(t) — Vi(t), and the overline denotes the time 
average. This correlation-based measure M c can be re- 
garded as a "statistical-mechanical" measure because it 
quantifies the average contribution of (microscopic) indi- 
vidual potentials to the (macroscopic) global potential. 
Hence, M c is in contrast to the conventional microscopic 
measure based on the cross-correlations between the indi- 
vidual potentials. As in the case of the whole population, 



the degree of inhibitory coherence in each subpopulation 
of the subthreshold and the suprathreshold neurons may 
be well quantified in terms of the coherence measures 
Mc and Mc pra ' based on the cross-correlations be- 
tween the global potentials (V su b and V SU pra) and the 
individual potentials, 



jLf(supra) 
c 

M {sub) 



1 



AT. 



supra . ^ 



]T Cf upra) (0), (15a) 



J-£c|^(o), 



(15b) 



C supra / \ 
i V) 



Ct ub (r) 



supra 

(t + r)Avi(t) 

^Avi pra (t) 



~~ 1 1 ••*) N supra): 



(16a) 



AV sub (t + t) Avi(t) 



^AVl b (t)^m 



(i = 1, ...,N sub ). 



(16b) 



By varying P SU pra , we measure the degree of inhibitory 
coherence in terms of the correlation-based measures M c , 
Mc Supra \ and Mc not only in the whole population, 
but also in the subpopulations of the subthreshold and 
the suprathreshold neurons, and the results are shown 
in Figs. BIbl)-|3Ib3). All of the coherence measures in- 
crease rapidly near the threshold Psupra > an d then they 
grow slowly. The values of these correlation-based mea- 
sures are very large in contrast to the spiking coher- 
ence measure M s . We also note that the degree of in- 
hibitory coherence in the subpopulation of subthresh- 
old neurons is higher than that in the subpopulation of 
suprathreshold neurons. This can be understood from 
the oscillating patterns of the global and the individ- 
ual potentials. The global potentials V SU pra and V su b 
exhibit small regular oscillations [see Figs. [3Ja2)-[2](a5) 
and Figs. [3Jc3)-[3Jc5)]. Like the case of the global po- 
tential, the individual subthreshold neurons exhibit only 
coherent subthreshold hoppings, in contrast to the case 
of suprathreshold neurons exhibiting both the coherent 
intermittent spikings and the coherent hoppings. Hence, 
the cross-correlations between V su b and the individual 
potentials of subthreshold neurons become higher than 



those between V s , 



and the individual potentials of 



suprathreshold neurons. 

In the above, we study the heterogeneity-induced in- 
hibitory coherence for a fixed value of D = 8 where 
Psupra — 0.16. By varying the noise intensity D we in- 
vestigate the effect of noise on the inhibitory coherence 
for J = 20. For P SU pra = 1, plots of the order param- 
eter O versus D are shown in Fig. 03a). The degree 
of inhibitory coherence decreases monotonically with in- 
creasing D from zero, and a transition to an incoherent 
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FIG. 5: (a) Plots of the order parameter O versus the noise 
intensity D for P 3up ra ~ 1, (bl)-(b5) raster plots of neural 
spikes and (cl)-(c5) global potentials Vg for various values 
of P S upra and D in a heterogeneous ensemble of N globally- 
coupled subthreshold and suprathreshold type-I ML neurons 
for J = 20 mS/cm 2 ; N = 10 3 in (bl)-(b5) and (cl)-(c5). The 
value of Ioc,i for each subthreshold (suprathreshold) type-I 
ML neuron is randomly chosen with a uniform probability in 
the range of {Idc~ A, Inc) [(Ihc, Ihc + A)] where P DC = 40 
/iA/cm 2 and A = 10 ^A/cm 2 . 



800 

200 
800 

200 
800 

200 



(a1) p supra =° 


(a2) 


(a3) p sup[a ^ 4 


(a4) 


(a5) p upra = 1 


_. i ■.. . i 




;*U!ti-£ : s '■:•'.« 
'.I 1 :" .."■i.vv "/ '. 


'J %\-*Mi 
[J- -i '■ \ \ i\: I 


| i"i j * i* li ( 1 1*" 










* i \ 
















• > j t -\ v 






iv. . " r 








... 



1200 1800 1200 1800 1200 1800 1200 1800 1200 1800 

t(ms) 

(b1 ) P, a „=0 (b2) P„„=0.2 (b3) P, upr ,=0-4 (b4) p ,„,„ =0 6 (b5) p ,„„= 1 



> -25 

E 







?WWVVWW\M 


viMMP/WW 










fwvwi/WWVW 































1200 1800 1200 1800 1200 1800 1200 1800 1200 1800 

t (ms) 




M «(1/M =1/M -1/N) 

syn.eff * Syn.eff syn ' 



FIG. 6: (al)-(a5) Raster plots of neural spikes and (bl)-(b5) 
global potentials Vg for various values of P sup ra and M syn in a 
heterogeneous ensemble of N (= 10 3 ) randomly-coupled sub- 
threshold and suprathreshold type-I ML neurons for J = 20 
mS/cm 2 and D = 8 /iA ■ ms 1 ' 2 /™ 2 . (c) Plots of the order 
parameter O versus the effective average number of synaptic 
inputs per neuron M syn]eff (l/M SJ(n>e // = l/M syn - 1/jV) 
for P suvr a = L The value of Inc,i for each subthresh- 
old (suprathreshold) type-I ML neuron is randomly chosen 
with a uniform probability in the range of (Ihc ~ A, Ihc) 
[(IdcJdc + A)] where P DC = 40 /iA/cm 2 and A = 10 
/xA/cm 2 . 



state occurs when passing a threshold D* (~ 28). Fig- 
ures EJbl)-[5jb5) and[5][cl)-[5][c5) show the raster plots of 
spikes and the global potentials Vq for various values of 

D and Psupra' 

For Psupra — 1, with increasing D the 
stripes in the raster plot become more smeared and the 
amplitude of Vg decreases. Eventually when passing the 
threshold D* , incoherent states appear (i.e., the raster 
plot consists of randomly scattered spikes and Vg ex- 
hibits a nearly stationary irregular oscillation). Hence, as 
D is increased the value of Pg Upra (threshold value of the 
fraction of suprathreshold neurons for the emergence of 
inhibitory coherence) increases; Psupra = 0-08 and 0.28 
for D = and 15, respectively. Thus, for D > D* no 
inhibitory coherence emerges, as shown in the case of 
D = 40. 

So far, we consider the globally-coupled case. How- 
ever, in a real brain each neuron is coupled to only a 
certain number of neurons which is much smaller than 
the total number of neurons. Due to the sparseness of 



the network architecture, the inhibitory coherence (seen 
in the globally-coupled case) is expected to be reduced 
or destroyed. It is often assumed in models that the cou- 
pling between neurons is random [lj [li| l4lTj43| . We 
briefly investigate the effect of sparse random connectiv- 
ity on the inhibitory coherence for J = 20 and D = 8 by 
varying the average number of synaptic inputs per neu- 
ron M syn in a heterogeneous ensemble of N randomly- 
coupled subthreshold and suprathreshold type-I ML neu- 
rons. Figs. [HIal)-[HIa5) and Etbl)-[HKb5) show the raster 
plots of spikes and the global potentials Vg for vari- 
ous values of M syn and P SU pra when N = 10 3 . For 
Psupra — 1> with decreasing M syn the stripes of spikes 
in the raster plot become more smeared and the am- 
plitude of Vg decreases. Eventually, incoherent states 
appear when passing a threshold M* yn (i.e., the raster 
plot is composed of randomly scattered spikes and Vg 
shows a nearly stationary irregular oscillation). Hence, 
as M syn is decreased from N — 1 (corresponding to the 
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globally-coupled case), a larger fraction of suprathresh- 
old neurons is necessary for the appearance of inhibitory 
coherence (e.g., see the cases of M syn =800 and 500). 
Thus, for M syn < M* yn inhibitory coherence disappears, 
as shown in the case of M syn = 100. As is well known, 
M syn [rather than P syn (i.e., the connection probability 
per neuron)] plays an appropriate sparseness parameter 
for the coherent transition because there exists a fixed 
threshold value M* yn for large N, independently of N 
2, 43]. (In contrast, the threshold value of P syn depends 
on N .) When P supra = 1, plots of the order parameter 
O versus the effective average number of synaptic inputs 
per neuron M syn ^ ef f (1/M syn , e ff = 1/M syn - l/N) are 
shown in Fig. Etc) for N = 100, 300, 500, 1000, and 3000, 
where the correction term (~ l/N) takes into account 
the finite network size effect @, [43[. Inhibitory coher- 
ence emerges when M syn e f f is larger than a threshold 
M* f , (~ 553); for the case of TV = 10 3 , M* yn ~ 356. 

Finally, we examine the universality of the 
heterogeneity-induced inhibitory coherence in a 
population of N globally-coupled subthreshold and 
suprathreshold quadratic integrate-and-fire (QIF) 
neurons |37| . The QIF neuron model is canonical in 
the sense that any type-I excitable systems can be 
transformed into the form of the QIF model by a 
continuous change of variables. Here, we use a QIF 
model with parameters derived from the Wang-Buzsaki 
conductance-based model [15]. These QIF neurons 
interact via the same inhibitory GABAergic synapses 
as those in the case of the ML neurons. The dynamics 
of the membrane potential V{ of each QIF neuron is 
governed by the following set of differential equations: 



c d ^ = A{v i -v*Y+i i -r 



(17) 

with the auxiliary after-spike resetting: if«j > vt (thresh- 
old potential for spiking), V{ <— v r (resetting potential), 
where t denotes the time (measured in units of ms), 
C = 0.9467 /iF/cm 2 , A = 0.012875 mS/cm 2 /mV, v* = 
-59.5462 mV, I* = 0.1601 ^A/cm 2 , v t = -26.3462 mV, 
and v r = —64.1462 mV. The ith neuron is stimulated 
by the applied current Ii and the white Gaussian noise 
£i whose intensity is controlled by the parameter D. The 
synaptic current injected into the ith neuron I syn ,i is 



given in Eq. (|2"c)) where the coupling strength is controlled 
by the parameter J, V syn {— —75 mV) is the synaptic re- 
versal potential, and the synaptic gate variable s obeys 
the 1st order kinetics of Eq. ([T};) along with Eq. ((2g) 
with 8 — 2 mV. For the single QIF neuron, a tran- 
sition from a resting to a spiking state occurs when Ii 
passes the rheobase I* via a saddle-node bifurcation on 
an invariant circle. Here a spread in the value of h is 
taken into consideration, and thus for each subthreshold 
(suprathreshold) neuron a value of U is randomly chosen 
with a uniform probability in the range of (I* — A,/*) 
[(/*,/* + A)], where the value of the spread parame- 
ter A is set as A = 0.1 /iA/cm 2 . By increasing the 
fraction of suprathreshold neurons P supr a , we investigate 
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FIG. 7: Order parameters, raster plots of neural spikes, and 
time series of global potentials in the heterogeneous ensem- 
ble of N globally-coupled inhibitory QIF neurons for J — 5 
mS/cm 2 and D = 0.1 ■ ms 1/2 /cm 2 ; N = 10 3 in (bl)-(b5). 
The value of Ii for each subthreshold (suprathreshold) QIF 
neuron is randomly chosen with a uniform probability in the 
range of (I* -A, I*) [(/*,/*+ A)] where P = 0.1601 fiA/cm 2 
and A = 0.1 ^A/cni 2 . Plots of the order parameter O versus 
the fraction of suprathreshold neurons Psupra in (a). Raster 
plots of neural spikes and time series of the global potential 
V G for Psupra = (bl) 0, (b2) 0.1, (b3) 0.3, (b4) 0.5 and (b5) 
1.0. 



the heterogeneity-induced inhibitory coherence for J = 5 
and D = 0.1. Figure [Tfa) shows plots of the order pa- 
rameter O versus P mmra - As P m 
value P*. 



supra • j supra passes a threshold 
,. upra {— 0.05), a transition from an incoherent 
to a coherent state occurs like the case of ML neurons 
[see Fig. (Hal)]. For an incoherent case of P SU pra = 0, 
the raster plot consists of randomly scattered sparse 
spikes and the global potential Vg exhibits a nearly 
stationary irregular oscillation, as shown in Fig. [TJbl) . 
However, when passing the threshold Ps Upra partially- 
occupied "stripes" appear in the raster plot along with 
regularly-oscillating small-amplitude global potential Vq 
[see Fig. Efb2) for P supra = 0.1]. As P sup ra is further 
increased, both the pacing degree of spikes and the am- 
plitude of Vg (representing the degree of inhibitory co- 
herence) increase, as shown in Figs. [7]Jb3)-[7][b5). 



IV. SUMMARY 

We have studied the heterogeneity-induced inhibitory 
coherence by increasing the fraction of suprathreshold 
neurons P sup ra in an ensemble of globally-coupled sub- 
threshold and suprathreshold type-I ML neurons. For 



supra 



= no inhibitory coherence has been observed, 
which implies that subthreshold type-I neurons are dif- 
ficult to synchronize by synaptic inhibition. However, 
as Psupra passes a threshold value Ps* upra , a coherent 
transition occurs in the subpopulation of suprathresh- 
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old neurons, and these synchronized suprathreshold neu- 
rons play the role of coherent inhibitors for the emer- 
gence of inhibitory coherence in the whole heterogeneous 
population. Consequently, for P sup ra > Psupra a regular 
population rhythm with small amplitude appears in the 
ensemble- averaged global potential Vg- For this coher- 
ent case, both the coherent intermittent spiking and the 
coherent subthreshold (small-amplitude) hopping phases 
appear in the individual potentials of suprathreshold neu- 
rons. Thus, these sparsely synchronized suprathresh- 
old neurons exhibit coherent mixed-mode oscillations 
with two well-separated frequency scales (i.e., fast sub- 
threshold frequency and lower spiking frequency). By 
virtue of their coherent inhibition, sparsely synchronized 
suprathreshold neurons suppress noisy activities of sub- 
threshold neurons. Hence, only the coherent fast sub- 
threshold hopping phase appears in the individual po- 
tentials of subthreshold neurons. The inhibitory coher- 
ence has been characterized in terms of the "statistical- 
mechanical" coherence measures based on spikes and cor- 
relations. The degree of inhibitory coherence was thus 
found to increase as P SU pra is increased for P sup ra > 
Psupra- The effect of sparse random synaptic connec- 
tivity on the inhibitory coherence has also been in- 
vestigated, and emergence of inhibitory coherence has 



thus been found to persist only if M syn is larger than 
a threshold value M* yn . Finally, the universality of 
heterogeneity-induced inhibitory coherence has been con- 
firmed in a population of canonical type-I QIF neurons. 
This kind of heterogeneity-induced weak inhibitory co- 
herence which emerges from sparsely synchronized os- 
cillations of suprathreshold neurons might be associated 
with cortical rhythms with irregular and sparse neural 
firings which contribute to cognitive functions such as 
information integration, working memory, and selective 
attention. Impaired sparse synchronization is believed to 
be related to mental disorders (e.g., schizophrenia and 
autism) @. 
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